NASA Technical Memorandum 105 762 


//760 3 



Computational Simulation of Matrix Micro-Slip 
Bands in SiC/Ti-15 Composite 


S.K. Mital 
University of Toledo 
Toledo, Ohio 

and 


H.-J. Lee, P.L.N. Murthy, and C.C. Chamis 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

1992 Winter Annual Meeting of the American Society of Mechanical Engineers 
Anaheim, California, November 8-13, 1992 


(NAS A- TM- 1057 6 2) COMPUTATIONAL N92-32466 

SIMULATION OF MATRIX MICRO-SLIP 
BANDS IN SiC/Ti-lS COMPOSITE 

(NASA) 9 p Unci as 




G3/24 0117603 



E-7179 


COMPUTATIONAL SIMULATION OF MATRIX MICROSLIP BANDS IN SiC/Ti-15 COMPOSITE 


S.K. Mital* 

University of Toledo 
Toledo, Ohio 43606 

and 

H.-J. Lee, P.L.N. Murthy, and C.C. Chamis 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

Computational simulation procedures are used to identify the 
key deformation mechanisms for a [0] g and a [90] g SiC/Ti-15 
metal matrix composite. The computational simulation procedures 
employed consist of a three-dimensional finite-element analysis and 
a micromechanics based computer code METCAN. The interphase 
properties used in the analysis have been calibrated using the 
METCAN computer code with the [90] 8 experimental stress-strain 
curve. Results of the simulation show that although shear stresses 
are sufficiently high to cause the formation of some slip bands in 
the matrix concentrated mostly near the fibers, the nonlinearity in 
the composite stress-strain curve in the case of the [90] g composite 
is dominated by interfacial damage, such as microcracks and 
debonding rather than microplasticity. The stress-strain curve for 
the [0] g composite is largely controlled by the fibers and shows 
only slight nonlinearity at higher strain levels that could be the 
result of matrix microplasticity. 

INTRODUCTION 

High temperature metal matrix composites (HTMMQ are 
potential structural materials for advanced propulsion systems 
demanding high operational temperatures (400 to 1100 °C). High 
stiffness to weight ratios, tailorable properties, dimensional stabil- 
ity, and hygral (moisture) resistance make these materials especially 
attractive for use in aerothermal structural applications. However, 
before HTMMC can be used in these critical applications, a thorough 
identification and understanding of their deformation behavior is 
required. 

Recently, matrix plasticity in the form of slip bands in SIC/ 
Ti-15 has been observed experimentally by Lerch (1991) using a 
post-heat treatment and etching method. These slip bands were 
observed in specimens strained to the nonlinear portions of their 
stress-strain curves but prior to failure. Results for [0] g and [90] g 


specimens with the fiber volume ratio of 0.34 show that slip bands 
are concentrated mostly around the Fibers and some interfacial 
debonding is observed for the [90] g specimen. However, no definite 
conclusions are made regarding the hierarchy of deformation. 

Similar studies conducted on SiC/Ti-15 by Newaz and 
Majumdar (1991) led them to conclude that plasticity is the dominant 
mechanism for the [0]g specimen. On the other hand, the [90] g 
specimen is controlled by both damage and plasticity. Newaz and 
Majumdar also suggest that the Poisson’s ratio can provide additional 
information regarding the amount of plasticity and damage that 
occurs in SiC/Ti-15. They suggest that an increase in major 
Poisson’s ratio is indicative of plastic deformation while a de- 
crease corresponds to the dominance of damage in the composite. 

The objective of this paper is to computationally simulate 
the stress-strain behavior of both [0] g and [90] g SiC/Ti-15 com- 
posites with a fiber volume ratio of 0,34 using METCAN and 
three-dimensional finite-element analysis (FEA). The resulting 
microstresses from both analyses are evaluated for the possibility 
of formation of matrix microslip bands and the key deformation 
mechanisms are identified. The computational simulation proce- 
dures consisting of a three-dimensional finite-element analysis and 
a micromechanics based computer code METCAN, are described 
briefly in the following sections. 


METCAN COMPUTER CODE 

METCAN stands for Metal Matrix Composite Analyzer and 
is a computer code developed at NASA Lewis Research Center to 
Simulate the high temperature nonlinear behavior of metal matrix 
composites. METCAN treats material nonlinearity at the constitu- 
ent (fiber, matrix, and interphase) level, where the behavior of the 
material is modeled using the dependence of a constituent’s prop- 
erties in a time- temperature- stress "material behavior space." The 
composite properties are synthesized from the constituent instanta- 
neous properties by virtue of composite micromechanics, compos- 
ite macromechanics, and laminate theory. Factors which affect the 
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behavior of the composite properties include the fabrication pro- 
cess variables, in-situ fiber and matrix properties, and/or the prop- 
erties of an interphase between the fiber and matrix. A unique 
aspect of the METCAN code is an integrated cyclic arrangement 
which defines the computational effort for each load increment as 
sh^wn in Fig. 1. Another feature is the use of a "multifactor 
interaction relationship" (MFIR) to represent the various nonlinear- 
ities and their interactions in the constituents. The multifactor 
interaction relationship (MFIR) and the reasons for its selection 
are shown in Fig, 2. More information regarding METCAN and the 
efforts to verify the code can be found in Chamis et al. (1990) and 
Lee et al. (1991). 


FINITE ELEMENT MODEL 

The finite element model used in this simulation consists of 
a group of nine fibers, all unidirectional, in a three- by-three unit 
cell array ("nine-cell model") as shown in Fig. 3. Since both the 
composites examined are unidirectional and to keep the time 
required in the analysis to a reasonable limit, only three plies (or 
three-by-three unit cell array) have been simulated in the finite- 
element analysis. The fiber diameter is taken as 150 pm and the 
interphase thickness is assumed to be approximately 3 percent of 
the fiber diameter with a Fiber volume ratio of 0.34. There are 16 
elements ("bays") along the length of the fiber. Each unit cell, as 
shown in Fig. 3, consists of 40 hexahedron (six-sided) solid ele- 
ments and eight pentahedron (five- sided) solid elements for a total 


of 6912 elements and 6953 nodes in the model. This corresponds 
to approximately 20 000 degrees of freedom in the model, which 
provides sufficient detail of the stress distribution in the composite 
for the purposes of this study. A general purpose finite-element 
analysis computer cock MSC/NASTRAN (1985) has been used to 
cany out the analysis. For longitudinal loading, uniform boundary 
displacements are specified in the 1-1 (along the fiber) direction. 
For transverse loading, uniform boundary displacements are speci- 
fied in the 2-2 (transverse to the fiber) direction. The resulting 
nodal forces to sustain the applied displacements are computed by 
the finite-element analysis and the resulting longitudinal or trans- 
verse stress is obtained by dividing the total force by the appropri- 
ate cross-sectional area. 


COMPUTATIONAL SIMULATION PROCEDURES 

The computational simulation procedures used in this study 
begin with the METCAN analysis to determine some of the con- 
stituents material properties for the finite-element analysis. Al- 
though bulk properties of the fiber and matrix are available (Lerch, 
1989) and listed in Table I, additional properties for the interphase 
are required to conduct an accurate finite-element analysis. How- 
ever, since experimental values for the interphase properties are 
unavailable, a calibration is conducted using METCAN with the 
[90] g experimental stress-strain curve to determine the interphase 
modulus and strength. Also, previous METCAN verification con- 
ducted on SiC/Ti-15 (Chamis et al., 1990) has shown that the 
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Rationale: 

• Gradual effects during most range, rapidly degrading near final stages 

• Representative of the in situ behavior for fiber, matrix, interphase, coating 

• Introduction of primitive variables (PV) 

• Consistent in situ representation of all constituent properties in terms of PV 

• Room-temperature values for reference properties 

• Continuous interphase growth 

• Simultaneous interaction of all primitive variables 

• Adaptability to new materials 

• Amenable to verification inclusive of all properties 

• Readily adaptable to incremental computational simulation 


■Matrix 



Notations: 

P - property; T - temperature; S - strength; R - metallurgical reaction; N - number of cycles: 
t - time; over dot - rate; subscripts: 0 - reference; F - final; M - mechanical; T - thermal 


Figure 2. — Assumed multifactor interaction relationship to represent various factors 
which influence in-situ constituent materials behavior. 



Unit cell 

691 2 Elements 
6953 Nodes 


Figure 3. — Schematic of model used for finite element analysis. SiC/Ti-1 5; fiber volume 
ratio, 0.34 
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Table 1: Constituent Material Properties. 



SiC Fiber 

Ti-15 matrix 

Interphase 
(3% thickness) 

Density, p, lb/in. 3 

0.110 

0.172 

0.172 

Stress free tempera- 
ture, T m , °F 

4870 

1400 

1400 

Elastic modulus, E, Msi 

62.0 

12.3 

17.0 

Poisson's ratio, v 

0.19 

0:32 

0.32 

Coefficient of thermal 
expansion, a, ppm/°F 

2.72 

4.72 

4.72 

Strength, ksi 

500 

60 

48 


matrix in the in-situ state will not realize the frill bulk matrix 
strength. Thus, the in-situ matrix strength is also calibrated from 
the [90] g stress-strain curve. 

The calibration procedure consists of three steps. In the 
first step, the in-situ matrix strength is varied until a match is 
obtained between the METCAN prediction for strength and the 
experimentally observed value for the [90] g laminate. The inter- 
phase modulus is determined in the second step by matching the 
initial slope of the METCAN predicted stress-strain curve with the 
experimental [90] g curve. In the third step, the interphase strength 
is calibrated by matching the first kink of the experimental [90] g 
curve. It should be noted that the calibration procedure is per- 
formed using only data from the [90] g case. Once the calibrated 
properties are obtained, they are kept at the same values through- 
out the study. 

After the calibration procedure is completed, three cases are 
examined using both METCAN and FEA. The first case simulates 
processing alone. The second and third cases simulate processing 
followed by a longitudinal loading on a [0] g and a [90] g compos- 
ite, respectively. Results from each case will be discussed in the 
following section. 

Processing Temperature: 

1292 to 70 °F 

Interphase Residual Stresses: 

Radial Hoop 

stress, stress, 

ksi ksi 

METCAN b -13 35 

FEA a -18 41 

a, b represent the location where 
the stresses are evaluated 


RESULTS AND DISCUSSION 
Processing 

Processing is simulated as a cool down from the processing 
temperature (700 °C/1292 °F) to room temperature (2\ °CP0 °F) 
by both METCAN and FEA. The METCAN analysis automatically 
takes into account the nonlinear variation of material properties 
with temperature and stress via the multifactor interaction relation- 
ship mentioned before. The nonlinear material behavior is incor- 
porated into FEA by using material properties generated from the 
multifactor interaction relationship. Ffer simplicity, the FEA simu- 
lation of processing has been carried out in one step using material 
properties computed at 360.5 °C/681 °F (average of 700 and 21 °C) 
by the multifactor interaction relationship using only a temperature 
dependence. Thus while the METCAN analysis takes into account 
the variation o f mat erial properties with temperature and stress, the 
constituent material properties depend only on temperature in the 
finite-element analysis. 

The residual stresses_that develop in the interphase as a result 
of processing at approximately the 3 o’clock position (point b) are 
shown in Fig. 4. FoFboth raSial and hoop stresses, METCAN 



METCAN unit cell 


Figure 4.— Results for processing. 
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predictions are conservative of FEA. This is expected since 
METCAN stresses represent average values over a whole region, 
while FEA determines stress values at a particular point. The 
large hoop stresses predicted by both analyses (35 to 41 ksi) due 
to processing alone are sufficiently high to cause the formation of 
radial microcracks in the interphase. Such radial microcracks have 
been observed experimentally by MacKay (1990). 

r901 g Composite 

The transverse stress-strain curves obtained from METCAN 
and FEA, along with the experimental curve are shown in Fig. 5. 
The curves are in good agreement. The first kink in the stress- 
strain curve occurs at an applied stress of approximately 42 ksi as 
interfacial debonding initiates. This applied stress level corre- 
sponds to a maximum shear stress of approximately 25 ksi. Al- 
though this shear stress is large enough to account for the local 
slip bands observed experimentally, examining the unloading curve 
in Fig. 6 shows that the nonlinear behavior is dominated by dam- 
age. Since there is good agreement in the loading part of the stress- 
strain curve (Fig. 5), only an average loading curve has been 
plotted in Fig. 6. The experimental modulus during unloading, as 
shown in Fig. 6, is significantly reduced in comparison to the 
loading modulus, which results in a very small permanent set 
(0.0004 total strain). The presence of a permanent set indicates 
that some plasticity is occurring, but the almost negligible value 
indicates that a damage mechanism plays a more dominant role. 
Since experimental observations have indicated both interfacial 
debonding (Lerch, 1991; Lerch and Saltsman, 1991) and interfacial 
radial cracks (MacKay, 1990) for this orientation, interfacial 
debonding is the critical deformation mechanism. 

This conclusion is reinforced by the predictions of major 
Poisson’s ratio by METCAN and FEA shown in Fig. 7. METCAN 
values are shown by the continuous line, while only two points are 
calculated from FEA. The decrease in Poisson’s ratio with increas- 
ing strain shown in the figure agrees well with the results observed 
by Newaz and Majumdar (1991), indicating that the primary defor- 
mation mechanism is damage. 


TOlg Composite 

The longitudinal stress-strain curves at room temperature 
obtained from METCAN and FEA are shown in Fig. 8, along with 
the experimental curve. The FEA curve is slightly on the higher 
(stiffer) side. This may be due to the assumption of only a tem- 
perature dependence of the constituent properties for FEA, while 
METCAN assumes a dependence on temperature and stress. The 
stress-strain curve is mainly linear with a slight nonlinearity at the 
higher strain levels. The maximum shear stress obtained from both 
FEA and METCAN is approximately 30 ksi. Once again, although 
this shear stress is sufficient to cause the formation of matrix slip 
bands observed experimentally, the mostly linear (elastic) behavior 
of the [0]g stress-strain curve indicates that plasticity is unlikely to 
be the dominant deformation mechanism. In order for any slip to 
take place, some type of damage in the form of either radial micro- 
cracks or debonding must exist in the interphase. In other words, 
the presence of such micro slip bands is indicative of interfacial 
damage, which could arise during fabrication. Previous work by 
the first author on composite microfracture (Mital et al., 1990) has 
shown that even an existing interfacial debond will not propagate 
if the loading is applied along the fiber (longitudinal) direction 
without a fiber fracture. Therefore, local interphase shear failure 
appears to be benign under this loading case and the composite 
behavior is controlled by the fibers. Thus, even with the presence 
of experimentally observed radial cracks in the interphase (reaction 
zone cracks) due to processing and the METCAN predicted failure 
of the interphase at a total strain of 0.006, only a slight nonlinear- 
ity is introduced in the stress-strain behavior after the interphase 
failure. This indicates that the observed matrix plasticity follows 
reaction zone cracks, and effects the stress-strain curve only at 
higher strain levels. This conclusion is further supported by the 
prediction of almost constant major Poisson’s ratio for this case as 
shown in Fig. 9. 
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Figure 7. — Variation of v 21 vs. strain at room temperature. 
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Figure 8— Room-temperature longitudinal [0] 8 stress-strain curve. 
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Figure 9.— Variation of v 21 vs. strain at room temperature. 


SUMMARY 

A three-dimensional finite-element analysis and the METCAN 
computer code were used to computationally simulate the stress- 
strain behavior of a [0] g and [90] g SiC/Ti-15 composite. Inter- 
phase properties and the in-situ matrix strength were determined 
by calibrating the METCAN using the experimental results for the 
[90] g composite. Results indicate that although matrix microslip 
bands are likely to occur concentrated mostly near the fibers, the 
nonlinearity in the composite stress-strain curve is dominated by 
damage either in the form of microcracks or debonding. In the 
case of the [90] g composite, the unloading modulus is much less 
than the loading modulus and the permanent set is extremely small, 
indicating the dominance of damage rather than microplasticity as 
the key deformation mechanism. In the case of the [0] 8 compos- 
ite, the stress-strain behavior shows only a slight nonlinearity at 
higher strain levels that could be the result of matrix microplastic- 
ity. The behavior of [0] g composite is controlled by the fibers in 
general and hence, any local shear failures are benign in this 
loading case. 
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